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The geometrical phase of a time-dependent singular oscillator is obtained 
in the framework of Gaussian wave packet. It is shown by a simple geo- 
metrical approach that the geometrical phase is connected to the classical 
nonadiabatic Hannay angle of the generalized Harmonic oscillator. 
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Explicitly time-dependent problems present special difficulties in classi- 
cal and quantum mechanics. However, they deserve detailed study because 
very interesting properties emerge when, even for simple linear systems, some 
parameters are allowed to vary with time. For instance, particular recent in- 
terest has been devoted to systems in which evolution originates geometric 
contributions [1-6]. One of these, the generalized harmonic oscillator has 
invoked much attention to study the nonadiabatic geometric phase for vari- 
ous quantum states, such as Gaussian, number, squeezed or coherent states, 
which can be found exactly [7-10]. Recently, the geometric phase for a cyclic 
wave packet solution of the generalized harmonic oscillator and its relation 



time-dependent Heller Gaussian wave packet form [11] 

®(x,t) = exp(ft _1 [-a(x - q) 2 + ip(x — q) + k]) (1) 

centred around the classical guiding trajectory (q,p), and proceed to derive 
equations of motion for the complex or real parameters (a(t), q(t), pit), and 
k{t)) which serve to specify a complete quantum wave packet. 

On the other hand, the number of exactly solvable quantum time-dependent 
problems is very restricted, one of the rare examples admitting exact solu- 
tions of the Schrodinger equation and have been studied intensively lately 
[12-19] is the quantum time-dependent generalized singular oscilator 

H = \ [Zit)p 2 x + Y(t)(p x x + xp x ) + X(t)x 2 + ^f-] (2) 

where x and p x = —ihd/dx are the quantum operators, X(t), Y(t), and Z(t) 
are an arbitrary function of time, and / is an arbitrary constant which could 
be zero. A distinguished role of the Hamiltonian (2) is explained by the 
fact that, in a sense, it belongs to a boundary between linear and nonlinear 
problems of classical and quantum mechanics. For this reason, it was used 
in many applications in different areas of physics. For example, it served as 
an initial point in constructing interesting exactly solvable models of inter- 
acting N-body systems [12-13]. It was also used for modeling diatomic and 
polyatomic molecules [14]. It can have some relation to the problem of the 
relative motion of ions in electromagnetic traps [19]. 

The aim of this letter is to explore Gaussian wave function dynamics 
for the Hamiltonian (2) with nonadiabatic time dependence, and formulate 



quantum and classical mechanics. For that purpose we intoduce a class of 
wave function of the form 

*,(*, t) = M-VmWTI7Z) exp | 1 (iM+ipq) [(^) 2 +2(^)] +k) } 

(3) 

given as products of squeezed Gaussian wave packet of the type (1) and a 



function x of order (1/2— sj (l/h) 2 + 1/4). Inserting Eq. (3) into Schrodinger 
equation 

ih-^ = ff*,, (4) 

and then compare between the coefficients of various powers of (x — q). This 
lead to 

(x - q) 2 : i$ = 2Z(3 2 - 2iY(3 + — , (5) 

Zi 

where (3 = -ip/(2q) - l/(2q 2 ), 



I Zl 2 
(x-q) 1 : -2i([3-—)(q-Zp-Yq) + (p + Xq + Yp-—)=0, (6) 
iq q 



[x — q)°: 

pq+tk = -[Zp 2 +2Ypq+Xq 2 + ^} ^L+2ft(l/2- ^JW + l/4)(Z0-i-). 

Z C[ C[ Z 

(7) 

The (x — q) 2 condition determines (3 by a nonlinear equation of the 
Riccati form, which can be transformed to a linear system by introducing a 
two dimensional vector v T = (Q,P) and 



where Q and P may be complex. In order that (5 satisfies (5), it is sufficient 
that v obey the Hamilton's equation 

-*)(?)-*»■ < 9 » 

The (x — g) 1 condition makes sense if 

q = Zp + Yq 

Zl 2 

p = -Xq-Yp+^ r , (10) 
Q 

and determines the complex guiding trajectory associated to the classical 
Hamiltonian 

H(q, p,t) = l [Z(t)p 2 + Y(t)(pq + qp) + X{t)q 2 + ^f-] (11). 
z q 

The (x — q)° condition determines the time dependent global phase and 
normalization included in k which can be rearanged in the form 

k(t)-k(0) = i f (L(t') + -2h(l- ^{l/hy + 1/4) (ZP - iY/2) ) dt', 

Jo ^ 1 ' 

(12) 

where L(t) = p(q,q)q — H[q,p(q, q),t\. Examining the three terms in the 
expression (12) for k, we see that the two first give i{(pq) — p(0)g(0)}/2. The 
remaining terms in k/h, namely -2i(l - ^{l/K) 2 + 1/4) J*(zp-iY/2)dt' is 
recognized as the angle (or phase) accumulated in the nonadiabatic evolution, 
and contain a dynamical part 9d(t) and a geometrical one 0jf(t) (Hannay's 
angle) . 

In the literature, 9n(t) is usually defined in relation with the intro- 



geometrical approach ( hence no calculatory) may be formulated to deduce 
the decomposition of the total phase factor. Let us consider the classical 
equation (9). The main property of this evolution is that is linear and area 
preserving. This implies that any initial conditions (Q(0),P(0)) being at 
t = on a centered ellipse £(0) in phase space evolve at time t on a similar 
ellipse £(t) of the same area. A little thought show that, more precisely, two 
points M , N on £(0) whose parameters differ by A(p evolve in points M u N t 
on £(t) with the same difference of parameters so that the area 
remains constant. The reason is that the standard parameter (p((p G [0, 2ir]) 
which parametrizes a point on an ellipse is such that it is proportional to 
the area swept by the vector v = OWL. The natural origins of these familly 
(homothetical centered) ellipses are the points associated with (p = 0, and 
the action angle coordinates (/, 9) of a point M in phase space, with respect 
to the familly associated with £{t), are respectively the area (divided by 2n) 
and the usual angular variable defined on this ellipse. There exist a natu- 
ral transport (with respect to the symplectic structure in phase space) for a 
familly of (homothetical centered) ellipses. More precisely, let M t on £{t) a 
point of coordinates (/, 0(t)) and M t +dt its transported on £(t + dt). Clearly, 
this transport must preserve areas, but this is not sufficient since it remains 
to precise how one point (for example the origin) is transported: for this, 
one simply requires that the area which is swept by the vectors v(ip) (on an 
ellipse) during the transport has, when averaged over <p, a mean value equal 
to zero. Therefore the transport is defined by 



This transport associates to M t of coordinates 0(t), the point Mt+dt of coor- 
dinates 0{t) +d0H(t) (0H(t)'- Hannay's angle). 

The angular coordinate on S(t + dt) of evolved point M t +dt of M t is 
0{t) + d0(t) . The difference d0(t) - d0 H (t) = d0 D (t) does not depend on 
the chosen point M t on £{t) . The quantity Id0o = I(d0 — d0n) is the area 



O Mt+dt /\ Mt+dt Mt+dt and can be written as the difference OM t /\M t M t +dt 



OM t A M t M t +dt of the area swept by the vector OM t during its evolution 
and of the one swept by OM t in the geometrical transport. Averaging these 
two areas on the ellipse £(£), one can see that Id0o can be interpreted as the 
mean value of the swept area during the motion. This justifies the appelation 
of dynamical angle for #d(£). 

To translate analyticaly the previous geometrical remarks, let £ be a 
complex two dimensional vector, it is known (for instance from optics) that 



one can describe (homothetical centered) ellipses as the set of vectors OM t — 
[E^e - ^] and those obtained by the transport are represented by the set of 



vectors OM t +dt = [E(t + dt)e~ lddH e~ %v }. Using the area preserving property 
(proportional to iE* A E) , one can easily verify that the transport defined 
by Eq. (13) is also written as E* A (dE - id0nE) = or 



3 E AE MA\ 

= — =5- • 14 

i(E* A E) 



Obviously Eq. (14) gives new expression for the nonadiabatic Hannay angle 
of the generalized harmonic oscillator. Within such formalism the above 



be looked for in the form 

v(t) = Ae-W> + rtE(t), (9(0) = 0), (15) 

with iE* A E conserved and A and ip fixed (A and ip are the conditions 
mesured with respect to the familly E(0)). The angular drift of the origins 
points of E(0) (measured with respect to E(t)) is naturally decomposed in a 
geometrical part (Hannay's angle) and a dynamical one 

i(E*AE) i(E*AE) 

(This relation is obtained by inserting v(t) in Hamilton equation (9) and 
making the wedge product with v*(t)). 

If one wants to explicitly calculate 9h and 9, we must find vector E(t). 
Before hand, one can impose the condition —i(E* A E) =41 (which fix the 
ellipse area) and take the first componant of E to be real. The general form 
of E(t) is thus 

g » = (^(&)' (17) 

It leads to the relations 

0=-2(Zl3-,Y/2)-~Ln(QQ'). (19) 

From the form of the time-dependent phase factors included in k Eq. 
(12), it is clear that 



(1 - y/{l/h)* + l/4) [A9(t) + l -Ln{ 



Q(t)Q*(t) 

Q(0)Q*(0) 



}]■ 



(20) 



we see that the logarithm term goes "downstairs" as time-dependent nor- 
malization factor in the ^i(x,t). The remaining term in 7/(t) is recognized 
as the phase factors acquiered by the wave packet in its evolution. 

Then, we can reach a simple relation between the geometrical phase for 
the quantum singular oscillator and the nonadiabatic Hannay angle associ- 
ated to the generalized harmonic oscillator 



where the first part is independent of h and is equal to the Hannay angle, 
the second part depends on h and /. 

Now we turn our attention to the Eq. (3) which can be rewritten as a 
simple Gaussian wave packet 



which is also a solution to the Schrodinger equation. 

In conclusion, the squeezed Gaussian wave packet dynamics for the time- 
dependent singular oscillator has been obtained as a simple squeezed wave 
packet. The quantum phases are obtained explicitly and connected to the 
classical angle of the generalized harmonic oscillator. The classical version of 



7 ?(f) = (1 - y/(l/h)* + l/4)A6 H (t), 



(21) 






(22) 



for the nonadiabatic Hannay angle has been obtained by employing a geo- 
metrical approach of the evolution in phase space. In the literature, 9n(t) is 
usually defined in relation with the introduction of time dependent canonical 
transformations. However, the geometrical part has been found by asking 
wether there exists a natural transport (with respect to the symplectic struc- 
ture) for a family of (homothetical centered) ellipses. When the parameter 
I vanishes, we see that the Gaussian wave packet (22) corresponds to the 
evolution of the "ground" state of the time-dependent generalized harmonic 
oscillator and the geometrical phase is equal to one-half of the classical angle. 
This is just what was obtained in Ref. [7]. 
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